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Abstract 


Theoretical expressions are derived for the steady-state frequency response 
of a thermocouple wire. In particular, the effects of axial heat conduction are 
demonstrated for both a uniform thermocouple wire and a non-uniform wire 
with unequal material properties and wire diameters across the junction. For 
the case of a uniform wire, the amplitude ratio and phase angle compare 
favorably with the series solution of Scadron and Warshawsky (1952) except near 
the ends of the wire. For the case of a non-uniform wire, the amplitude ratio at 
low frequency 0 agrees with the results of Scadron and Warshawsky for a 
steady-state temperature distribution. Moreover, the frequency response for a 
non-uniform wire in the limit of infinite length l is shown to reduce to a 
simple expression that is analogous to the classic first order solution for a 
thermocouple wire with uniform properties. 




I. INTRODUCTION 


The evaluation of jet engine performance and fundamental studies of 
combustion phenomena depend on the measurement of turbulent fluctuating 
temperatures of the gas within the engine (Dils and Follansbee, 1976). 
Historically, these temperatures have been measured with thermocouples. The 
design of a thermocouple represents a compromise between accuracy, ruggedness 
and rapidity of response. 

For example, the measurement of fluctuating temperatures in the high 
speed exhaust of a gas turbine engine combustor is required to characterize the 
local gas density gradients or convective heat transfer (Fralick, 1985). Although 
thermocouples are suitable for the measurement of high frequency temperature 
fluctuations (< 1 KHz) in a flowing gas or liquid, the measured signal must be 
compensated since the frequency of the time dependent fluid temperature is 
normally much higher than the corner frequency of the thermocouple probe 
(Scadron and Warshawsky, 1952). Moreover, the amplitude and phase angle of 
the thermocouple response may be attenuated by axial heat conduction for 
rugged thermocouples of finite length (Elmore, et al.; 1983, 1986). 

In the present study, the theoretical steady-state frequency response of a 
thermocouple wire has been calculated to include the effects of axial heat 
conduction. Solutions are derived for both a uniform wire and a non-uniform 
wire with unequal properties across the junction (e.g., different wire diameter. 
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thermoconductivity, etc.). Solutions are presented in the form of the amplitude 
ratio and phase angle for both the uniform and non-uniform thermocouple 
wires. In both configurations, the results are presented for thermocouple wires 
of either finite length where the effects of axial heat conduction are 
demonstrated or wires of infinite length for which the neglect of axial 
conduction is justified. For the case of the uniform wire, the theoretical results 
are compared with the series solutions of Scadron and Warshawsky (1952). 
Moreover, for non-uniform wires the results are compared in the limit of low 
frequency with the steady-state temperature solution of Scadron and 
Warshawsky. 


H. THEORY 


The steady-state frequency response of a thermocouple wire will be 
developed in the following sections with the following assumptions: (a) the 
amplitude of the fluctuating fluid temperature is small relative to the mean 
absolute temperature (b) the thermocouple dimensions are small relative to the 
size of the turbulent eddies or enclosure dimensions (c) radial temperature 
gradients in a wire cross section can be neglected and (d) radiative heat transfer 
can be neglected relative to conduction and convection. 

1. Uniform Wire 

In this section the geometry of Fig. 1 is considered where the material 
properties of thermal conductivity k, specific heat c and wire density pare 
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assumed to be equal on both sides of the thermocouple junction. If the probe is 
subjected to a flowing fluid, whose temperature fluctuates about its mean, the 
expression for the local conservation of energy in the thermocouple wire 
becomes (Scadron and Warshawsky, 1952) 

HI*. = + —It -T w ) (1) 

dt dx 2 pcD { ‘ w) 


where a = k/pc is the thermal diffusivity of the wire, T g is the ambient fluid 
temperature, h is the convective heat transfer coefficient, D is the wire diameter 
and T w is the local wire temperature measured along the axis at a distance x from 
the centerline of Fig. 1. 

If the wire and fluid temperatures are measured relative to the support and 
mean fluid temperature To, the variation in the ambient fluid temperature can 
be written in the complex form 

T t {t) = T 0 + T I e i °“ (2) 

where co is the angular frequency of the ambient fluid. Since Eq. (1) is linear, we 
now seek a solution for the local wire temperature of the form (Hildebrand, 1976) 

T w = T 0 + T u (x)e ia> '. (3) 

Referencing all temperatures with respect to the steady-state temperature T 0 and 
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normalizing with respect to the amplitude of the fluctuating ambient fluid 
temperature Tf, one defines a local normalized steady-state frequency response 

f( x) for the thermocouple wire of the form 


7>r 0 

T f 


= T ( x)e--. 


(4) 


Substituting Eqs. (2) and (3) into Eq. (1), one obtains an ordinary differential 
equation of the form 


icoT 


a LL + coS\~T) 

dx 


(5) 


where T = T (a:) is the frequency response and the vector notation will be dropped 
for simplicity. Thus for the geometry of fig. 1, one seeks a solution to the non- 
homogeneous linear second order differential equation for the dependent 
variable T of the form 


yT" - G(co)T = -1 


( 6 ) 


subject to the boundary conditions 


T(l) = T(-l) = 0. 


( 7 ) 


Here, the parameters in Eqs. (5) and (6) are defined as 
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( 8 ) 


co. = 


Ah 

pcD' 


a 



G(o)) = 1+ i(— ) 
co. 


where is the natural frequency of the uniform wire of fig. 1. 

The general solution to Eq.(6) can be written in the form (Hildebrand, 1976) 

T = A sinh qx + B cosh qx + } - (9) 

G((0) 

where the constants A and B are complex, l/G(<o) represents the particular 
solution and the parameter 

Substituting x = ±l in Eq. (9), the boundary condition Eq. (7) yields values for the 
constants 


A = 0, 


B = 


-1 

G(w)cosh ql 


( 10 ) 


Hence, one obtains a steady-state frequency response for the simple uniform wire 
of Fig. 1 in the form (Forney, 1988) 


T = 


G(OJ) 


(1 - 


cosh qx 
cosh ql 


(ID 
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The steady-state frequency response at the junction x = 0 of fig. 1 now becomes 


T(0) = — (1 " sech< 7 /). (12)' 

G(o» 

The term 1 - sech (q 1 ) in Eq. (12) corresponds to the steady-state temperature 
distribution for a constant ambient fluid temperature Tf > T 0 derived by 
Scadron and Warshawsky (1952). It is also interesting to note that the first order 
response can be derived from Eq.(12) in the limit of large wire length / or 

limT(O) = — (13) 

G(co) 


where Eq. (13) represents the frequency response of a uniform thermocouple 
wire with no axial heat conduction. 

The steady-state frequency at the thermocouple junction x=0 is normally 
characterized graphically in the form 

T( 0) = |r(0)|e'* (14) 


where \T(0)\ is the amplitude ratio and is the phase angle. In the latter case, the 
phase angle in degrees is 


O = 57.3 tan -1 


Im T(0) ~ 
Rer(O) 


(15) 
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where Im[T(0)] and Re[T(0)] are the imaginary and real parts of T(0), respectively. 


2. Previous Work 

Scadron and Warshawsky (1952) in their pioneering work considered the 
frequency response of a uniform thermocouple wire as shown in fig. 1. 
Introducing transformations for the dependent variable T w in Eq. (1) as suggested 
by Carslaw (1945), Scadron and Warshawsky obtained expressions for the steady- 
state frequency response in the form of the amplitude ratio 


|r(0)| = 


Is , 2 




Si 




( 16 ) 


and phase angle 


O = tan _1 (— ^r) _ tar*"' ( — )• 
to S. co. 


(17) 


Here, the infinite series Si and S 2 that appear in Eqs. (16) and (17) are 
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(18a) 


s, = i - Z 

>«1 


<W-» 




(0 

CD. 


I 


and 


^2 - X c \ 

}• i 


d y -l 


*/ + 


(0 


(18b) 


where the constants that appear in the terms of Si and S 2 are 

Cj = ^[l-cosC/g)]sin^ ^~ X) j (19a) 


and 


*j = 1 + 


U*rr 

4/ 2 ‘ 


(19b) 


Properties of these solutions derived by Scadron and Warshawsky are discussed 
in Sec. HI. 

3. Non-Uniform Wire 


If the wire diameters or the physical properties such as thermoconductivity 
are different across the thermocouple junction as shown in fig-2, Eq.(6) must be 
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solved in the separate regions 1 and 2 on both sides of the junction. This 
solution is subject to the boundary conditions of continuous temperature and 
equal rates of conductive heat transfer at the junction x = 0 or 


T x (-l) = r 2 (/) = 0 (20a) 

r,(0) = T 2 ( 0) (20b) 




d7;(0) 

dx 




dT^ 0 ) 

dx 


(20c) 


In region 1, the general solution to Eq. (6) for the frequency response can be 
written in the form 


T x (x) = \ sinh q x (/ + x) + B x cosh q x (/ + x) + 


1 

G x {(0) ' 


( 21 ) 


Since T x (-l) = 0 from boundary condition Eq. (20a), one finds for the constant Bi 


0= B x + 


G x (ca) 


Substituting for Bi into Eq. (21), one obtains 
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r,U) = A sinh q x (l + x) + ^-^[l-cosh< 7 ,(/ + x)] 


( 22 ) 


where the parameters in Eq. (22) reflect the material properties and wire diameter 
of region 1 or 


P\ C \D\ 



C,(co) = l + i 



Similarly in region 2, the frequency response can be shown to be of the form 


T 2 (x) = A 2 sinhq 2 (l-x) + —^—[l-coshq 2 (l-x)} 

U 2 (Ci)) 


(23) 


where the constants Ai and A 2 in Eqs. (22) and (23) above must be determined by 
the simultaneous solution to the boundary conditions Eqs. (20b) and (20c). 

Substituting Eqs. (22) and (23) into boundary condition (20b), one obtains a 
linear equation for the constants of integration Ai, A 2 or 

A sinh q x l - A sinh q 2 l = R x . (24) 

Similarly, substituting Eqs. (22) and (23) into boundary condition (20c), one 
obtains a second linear equation for Ai , A 2 or 

A Qcosh qj + ^coshqj = A . (25) 
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Here, the complex constant Q is defined as 


Q = 


ja A ffi 


( 26 ) 


and the complex constants Ri , R 2 are equal to 


A 


(-L - -L) + i- cosh - rq-cosh q 2 l 

G 2 G, Gj g 2 


(27) 


and 


A 


— sinh <7,/+— sinh q 2 l. 
Gj G 2 


(28) 


Solving Eqs. (24) and (25) for Ai and A 2 , one obtains the determinate system 


A « 


R { - sinh <y 2 / 
R 2 cosh q 2 l 

DET 


(29) 


and 


A = 


sinh <7,/ /?, 

Q cosh <7,/ /? 2 

DET 


(30) 
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where the determinate in the denominator is equal to 


DET = sinh q x l cosh q 2 l + Q cosh q x l sinh q 2 l. 


( 31 ) 


Thus, the constants Ai and A 2 in Eqs. (22) and (23) become 


sinh j^tanh ^ +Q tanh q 2 l^+ tanh ^ + tanh q 2 l^j 

** = ' sinh qj + Q cosh q x l tanh q 2 l 


(32) 


and 


^ = 


s^^tanh qj - lanh^ ^) + 2!|^(tanh 2' tanh q,l 

sinh q 2 l + Q' cosh q 2 l tanh qj 


(33) 


where the complex constant in A 2 is 


Q' = 1/(2- 

Finally, the steady-state frequency response at the thermocouple junction x 
= 0 can be determined by substituting the value for Ai above or Eq. (32) into the 
general solution for Ti(0) given by Eq. (22). One obtains after simplification 

g^Mfl-coth q,l tanh qj tanh f) 

Tm , . G. V 2 l ^ (34) 

i i q coth q x l + coth q 2 l 
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where it can be shown that Ti(0) = T 2 (0) at the thermocouple junction. 

It is interesting to note that one can determine a frequency response for a 
non-uniform wire of infinite length from Eq. (34) in a manner similar to the 
derivation of Eq. (13). Taking the limit of infinite wire length in Ti(0) above , 
one obtains the simple expression 


limr.(O) = 

/- M- 


Gj 

(2 + 1 


(35) 


Equation (35) represents the frequency response for a non-umform wire which is 
analogous to the well known first order response given by Eq. (13) for a uniform 
wire of infinite length. 


m. RESULTS 

The amplitude ratio and phase angle of the thermocouple frequency 
response were plotted graphically for the case of a uniform wire with the 
geometry shown in fig. 1. In this case, average properties of a type B or 
Pt 16% Rh-Pt / 30% Rh were used since the material properties were nearly equal 
across the thermocouple junction. The wire dimensions, properties and gas 
conditions are listed in table KTouloukian et al., 1970). 


The amplitude ratio and phase angle were also plotted for a composite or 
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non-uniform thermocouple wire with the geometry shown in fig. 2. In this case, 
different wire diameters and material properties were used including a type B 
thermocouple described in Table 1 in addition to a type T or copper - constantan 
described in table 2 (Touloukian et al., 1970). 

The form of the convective heat transfer coefficient h that appears in the 
computation of the natural frequency defined in Eq. (8) was determined from 
the expression (Scadron and Warshawsky 1952) 

Nu= .431 Re^ (36) 

where Nu(=hD/k f ) [ s the Nusselt number, k f is the thermal conductivity of the 
ambient fluid and Re (= v £)/ v f ) j s the Reynolds number of the thermocouple 
wire. Here, v and are the fluid velocity and kinematic viscosity, respectively. 

It should be noted that the convective heat transfer coefficient h a D' n and the 
natural frequency of a thermocouple wire for given material properties 

clD~ v \ 

1. Uniform Wire 

The amplitude ratio |7X0)| at the wire junction for the steady-state frequency 
response Eq. (12) is computed in fig. 3. This assumes a type B thermocouple wire 
with the dimensions listed in table 1. In this case, the average material 
properties listed in table 1 were used since the theory assumes that the 
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thermocouple is uniform across the junction. In fig. 3 the amplitude ratio |T(0)| 

-> .54 for low frequencies co -> 0. This lower limit corresponds to the attenuation 
expected for a thermocouple wire exposed to a steady ambient gas temperature Tf 
> To as derived by Scadron and Warshawsky (1952). 

The amplitude ratio Eq. (16) predicted by Scadron and Warshawsky (1952) is 
also graphed in fig. 3. It was found that the series Eqs. (18a) and (18b) converged 
slowly such that the final result was in error by ~ 11% at low frequencies after 
retaining the first 10 terms. The first 50 terms of the series expressions for Si and 
S 2 reduced the error to 2.3% while retaining 100 terms as was done in fig. 3 
limited the error to 1%. The final result is identical, however, to the amplitude 
ratio predicted by Eq. (12) and derived in the present study. 

The amplitude ratio for a type B thermocouple wire of infinite length from 
the first order response Eq. (13) is also indicated in fig. 3. As expected, the first 
order response or Eq. (13) asymptotically approaches the amplitude ratio 
predicted for a wire of finite length in the limit of large frequencies oj/o) m » 1 
where the time scale of the ambient temperature fluctuation is too small for axial 
heat conduction to occur. 

The phase angle o for all three theories is computed in fig. 4. The phase 
angle varies over the range 0 ^ <X> > - n/2 and approaches the lower limit of 
- n/2 for large frequencies w / w n » 1. It is apparent that the effect of axial 
heat conduction is to reduce the magnitude of the phase angle. 
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The spacial variation of the amplitude ratio |F(*)| derived from Eq. (11) is 

graphed in fig. 5. The present theory approaches the correct limit of |T'(±/)| = 0. 
However, the theory of Scadron and Warshawsky oscillates near the end of the 
wire x->±L In the latter case the amplitude ratio equals the value 

1 / *^1 + (<y/co.) for x =±l since S] = 1 and S 2 = 0 in Eqs. (18a) and (18b), 
respectively. An incorrect value for the amplitude ratio near the ends of the 
wire would not ordinarily present a problem unless one wished to adapt the 
theory to the common configuration of a supported thermocouple wire with 
heat transfer down the support legs (Fralich and Forney, 1990). 


2. Non-Uniform Wire 


The amplitude ratio fc(0)| from the steady-state frequency response Eq. (34) 
is graphed in fig. 6. . This plot assumes a type T thermocouple wire with the 
dimensions listed in table 2 . An important parameter that affects the relative 
magnitude of the terms in Eq. (34) and thus the relative importance of the 
thermocouple wire segments in regions 1 or 2 in the schematic of Fig. 2 is the 

complex dimensionless ratio Q = K D U\ / We note here that Q assumes 

different values in the limit of either low or high frequencies. Given that the 
Nusselt number is related to the wire Reynolds number by Eq. (36), one obtains 
for low frequencies 


Q{c* = 0) = Qo 






(37) 
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while for high frequencies 


Q(a) = °°) = (L = 


k 

1*2 


'it 

V^J 


(38) 


The amplitude ratio |7](0)|-> .228 for low frequencies a>-»0 in fig. 6. The 
lower limit corresponds to the attenuation predicted by Scadron and 
Warshawsky (1952) for a non-uniform thermocouple wire as shown in fig. 2 that 
is exposed to a steady-state ambient gas temperature T f > T 0 . The real 
expression for the temperature distribution developed by Scadron and 
Warshawsky for this case is 


IMzl o = 
Tf 


(39) 


where 


ly* [y~ 

w = sech cosh - J — \ 
r 2 2 1 


ly~ ly* 

1 tanh-^-coth— 

1 - m i 2 sinh ly~ csch ly* 


(40) 


Here, the parameters y + and y are defined as 



(41a) 


and 
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(41b) 


m, 2 


0,-1 

0) + i ' 


The amplitude ratio |71(0)| from Eq. (35) for a type T thermocouple wire with 
the same wire dimensions D 2 / Di * 2 as the solid curve in fig. 6 but of infinite 
length /-> 00 is also computed. As indicated, the amplitude ratio at low 
frequencies for a non-uniform thermocouple wire of infinite length 

approaches |r,(0)|->L This result is identical to the first order frequency response 
Eq. 13 for a uniform wire of infinite length as is shown in fig. 3. For the non- 
uniform wire, however, axial heat transfer does occur at the junction because of 
dissimilar material properties. At large frequencies o)fo) m » 1, however, the 
amplitude ratio derived from Eq. (35) for the wire of infinite length 
asymptotically approaches the value predicted by Eq. (34) for a wire of finite 
length as shown in fig. 6. The latter result is similar to the features of fig. 3 at 
large frequencies for the case of a uniform thermocouple wire. These latter 
features reflect the fact that the time scales of the ambient gas temperature 
fluctuations at high frequency are too small for axial heat conduction to occur. 

The phase angle <j> for both expressions Eqs. (34) and (35) is graphed in fig. 7. 
The phase angle varies over the range 0 ^ ^ — tt/ 2 and approaches the lower 

limit of - tf/2 for large frequencies . It is apparent that the effect of axial heat 
conduction is to reduce the magnitude of the phase angle. 

The amplitude ratio |T,(0)| derived from Eq. (34) for a type B thermocouple 
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wire is computed in fig. 8. Since the material properties in this case are nearly 
equal on both sides of the junction, the indicated difference in the amplitude 
ratio is the result of different wire diameters across the junction. These results 
are strongly influenced by the magnitude of the factor Q 0 defined by Eq. (37). If 
the parameter Q 0 ~ 1 as is the case for the upper curve in fig. 8, both regions 1 and 
2 in the schematic of fig. 2 are equally weighted in the expression for the 
frequency response given by Eq. (34). However, increasing the wire diameter D 2 
in region 2 means Qo is reduced in magnitude such that Qo < 1- Now the 
amplitude ratio is more strongly influenced by the properties of the larger wire 
in region 2. The phase angle is plotted in fig. 9 for the same wire configuration. 

The amplitude ratio ft Ml is also computed as a function of the axial 
distance down the type B thermocouple wire in fig. 10. This was accomplished 
by substituting the expression for the constant Aj given by Eq. (32) into the 
general solution for the frequency response Eq. (22). As discussed above, if the 
ratio Qo < 1 because of an increase in D 2 , a large increase in the rate of heat 
transfer occurs in region 2 down the larger wire. This additional axial heat 
conduction reduces the amplitude ratio at the thermocouple junction compared 
to the case of Qo - 1 for the wire geometry represented by the upper curve in fig. 

10 . 


The amplitude ratio derived from Eq. (34) for a type T thermocouple wire is 
graphed in fig. 11. Since the thermoconductivity of copper in region 1 is more 
than an order of magnitude larger than the thermoconductivity of constantan in 
region 2, the parameter Q 0 > 1 for Di/D 2 = 1. Increasing the diameter of region 2, 
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however, has a small effect on the amplitude ratio since Qo > 1/ which remains 
unchanged, and the properties of region 1 continue to dominate the rate of heat 
transfer. The phase angle for the same configuration is computed in fig. 12. 

The amplitude ratio fcWl is also computed as a function of the axial 
distance down the type T thermocouple wire in fig. 13. As discussed above, Qo > 
1 for both wire configurations. Thus, the large rate of axial heat transfer down 
the copper wire in region 1 gives nearly equal spatial temperature profiles for 
both wire configurations and this fixes the junction temperature. 
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IV. NOMENCLATURE 


A = 
B = 

c = 

q = 

dj = 
D = 

G = 

h = 

i = 

k = 

k f = 
/ 

mi2 = 
Nu = 

Q = 
Qo = 
Q- = 
q = 

Re = 
Si = 


constant of integration 
constant of integration 

material specific heat {j ~8 m ~ K 
term in infinite series Si , S2 
term in infinite series Si , S2 
thermocouple wire diameter (cm) 

l+ittu/co.) 

heat transfer coefficient [J -cm 1 - s 


unit imaginary number 


(=^i) 


material thermoconductivity { J cm 

gas thermoconductivity {j -cm - s 
length of thermocouple wire (cm) 

(a-i)/(a+i) 

Nusselt number (= hD/kf) 



{K!K)\DJDr 

( G/y ) 5 

Reynolds number (=vD/u f ) 
infinite sum 


21 



infinite sum 


time (s) 

frequency response 

amplitude of periodic gas temperature ( %) 
gas temperature 

steady-state gas temperature (°^) 

complex amplitude of periodic wire temperature 


local wire temperature (°^j 

gas velocity ( cm ~ s ' ) 

axial distance from center of wire (cm) 

Greek Symbols 

thermal diffusivity (cm 2 - s ') 
a/o) m (cm 2 ) 

Xtr, * y-m | ~" 1 

/vs - /« 

kinematic viscosity of gas ( cm * ~ 5 ') 

angular frequency (* ') 

natural frequency of wire ( = I pcD)(s l ) 

phase angle (degrees) 

normalized temperature at junction 

-3 

material density (gm - cm ) 
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VI. TABLES AND FIGURES 


Table 1 - Uniform Wire Conditions 
Dimensions (cm) Air Properties 

D = .0076 T 0 = 300 °K 


Pt - 6% Rh 
Pt - 30%Rh 
Average 


P = 1 atm 

v = 50 m/sec 

Properties of Type B 


/ ] 

ft 


[cm 3 -°Kj 

IX 

^ sec j 

2.73 

.238 

2.86 


.190 

2.8 

.214 


Table 2 - Properties of Non-Uniform Wire 
Dimensions (cm) Wire Location 


r 


Di(cm) 

D?(cm) 

EhZEh 

/(cm) 

Region 

Type B 

Type T 

.0076 

.0076 

1 

.1 

h 

Pt-6%Rh 

Copper 

.0076 

.0152 

2 

.1 

2 

Pt-30%Rh 

Constantan 


Properties of Type T 



f J ) 

[cm 2 ] 


pc\ 3 

a\ 


(^cm — K J 

L sec J 

Copper 

3.44 

1.16 

Constantan 

3.48 

.067 
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Fluid Flow 



Fluid Flow 



Fig. 2 Schematic of non-uniform wire 
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Type B (I//Y = 1.2) 
theory 

Scadron and Warshawsky (1952) 
first order (\/fy = oo) 


10 

“/“n 

Fig. 3 Amplitude ratio vs angular frequency for uniform wire. 
Solid line, Eq. (12); Circles, Eq. (13). 






x/l 

Fig. 5 Amplitude ratio vs axial distance for uniform wire. 
Dotted line, Eq. (11). 










Fig. 9 Phase angle vs angular frequency for non-uniform wire with similar 
material properties. Curves derived from Eq. (34). 
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Fig. 10 Amplitude ratio vs axial distance for non-uniform wire with similar 
material properties. Each curve is derived from Eqs. (22) and (32). 
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Fig. 13 Amplitude ratio vs axial distance for wire with dissimilar material 
properties. Each curve is derived from Eqs. (22 ) and (32). 




VII. FUTURE WORK 


The experimental apparatus to test various thermocouple configurations 
is fully constructed. Dr. Earl Meeks of the Georgia Tech Microelectronics 
Laboratory has been given a subcontract to construct the thermocouple probes. 
During the next 6 months, several thermocouples will be constructed and tested. 
These thermocuples will consist of either type K (chromel-alumel) or type T 
(copper-constantan) junctions and will have diameter ratios of either 2 or 5. 
Among these thermocouples will be several with support legs and some 
without. These data consisting of amplitude ratio and phase angle will be 
correlated with various theoretical expressions of the type presented in the 
present report. 
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